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The number of gluons in the hadron wave function is discrete, and their formation 
in the chain of small x evolution occurs over discrete rapidity intervals of Ay ~ l/a^ . 
We therefore consider the evolution as a discrete quantum process. We show that the 
discrete version of the mean-field Kovchegov evolution equation gives rise to strong 
fluctuations in the scattering amplitude, not present in the continuous equation. We 
find that if the linear evolution is as fast as predicted by the perturbative BFKL 
dynamics, the scattering amplitude at high energies exhibits a chaotic behavior. As 
a consequence, the properties of diffraction at high energies become universal. 

The wave function of an ultra-relativistic hadron is described by the Color Glass Con- 
densate (CGC) quasi-classical non-Abelian Weizsacker- Williams field ISL M . 
It emerges when the occupation number of the bremsstrahlung gluons emitted at a given 
impact parameter exceeds unity and eventually saturates at ~ l/os. It has been argued that 
in a big nucleus, such that asA^^^ ^ 1, and not very high energies the mean-field treatment 
is a reasonable approximation of the CGC evolution equations. In this approximation the 
interaction of a dipole with a big nucleus is decribed by so-called fan diagrams However 
in general the quantum fluctuations around the classical solution can strongly modify the 
scattering amplitudes. One way to understand the origin of these fluctuations is to consider 
the gluon emission as a discrete process It was argued in jsl that a color dipole wave 
function is dominated by rare dipole configurations which lead to strong fluctuations of the 
dipole (gluon) density and, as a consequence, to the fluctuations of the scattering amplitude. 
It was proposed to take such fluctuations into account by introducing new terms into the 
Kovchegov equation 0, Q . 

In this paper we take a different approach. An introduction of an infrared cutoff A on 
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the momentum of the emitted gluons amounts to imposing the boundary condition. This is 
equivalent to the quantization of the gluon modes in a box of size L ~ A^^, in which case 
the spectrum of the emitted gluons and their number become discrete. The formation of 
a gluon occurs over a rapidity interval of Ay ~ l/a^. Therefore, the evolution in rapidity 
can be considered as a discrete quantum process, where each subsequent step occurs when 
Ay as ~ 1. We will show that the discrete version of the mean-field Balitsky-Kovchegov 
evolution equation 3; ^3 gives rise to chaotic behavior of the scattering amplitude. We find 
that even the event-averaged scattering amplitude differs significantly from the continuous 
result. We will argue that this chaotic behavior is a general feature of discrete evolution 
when the growth of the scattering amplitude is sufficiently fast; this appears to be the case 
for the perturbative BFKL evolution. 
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In the conventional mean-field approximation Kovchegov equation jlll . 
for the scattering amplitude N{xq — Xi,b,y) of the color dipole of transverse size Xq — x^, 
at impact parameter b = (xg + Xi)/2 and rapidity y. In this paper we limit ourselves to the 
case of the most central collisions in which case the scattering amplitude can be evaluated 
at a fixed impact parameter. Hence omitting b the Kovchegov equation reads 

= ^ / d'x2^^iNix,„y) + iV(xo2,y) - Nix,„y) - N{x,,,y) N{x^,,y)) 
oy m J Xi2 

(1) 

It will be more convenient for our discussion to rewrite this equation in the momentum 
space. Defining 

N{k,y) = [ d'x\e'^-^N{x,y) (2) 



we can write (jH) as 

dN{k,y) 



x^ 



asX{l)N{k,y) - asN\k,y), (3) 



dy 

where x{l) is the familiar leading eigenvalue of the BFKL equation 

X(7) = 2^(1) - ^(1 - 7) - ^(7) (4) 
and the anomalous dimension operator is given by 

7W = l+33^. (5) 

Eqs. (H} and Q are written for the event averaged scattering amplitude. Let us, in the 
spirit of Ref. j7|, consider the evolution in a given event. The emission of one low x gluon is 
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FIG. 1: The scattering amplitude (scaled by — 1)) as a function of the number of evolution 
steps. The dashed curve indicates the continuous limit; the discrete evolution (solid line) leads to 
the chaotic behavior. 

of the order otsy ~ 1 contribution to the scattering amphtude. Therefore, we can enumerate 
the emitted gluons by discrete values n of the "time" parameter oisU- Thus the discrete 
version of Kovchegov equation takes the form 

Nn+i{Ly) = (1 + xmNn{k,y) - N^{k,y), (6) 

Let us expand x near the saddle point 70. In this paper we keep only the first term in 
the expansion for the sake of simplicity - this toy model suffices to illustrate our main idea. 
Introducing the rescaled amplitude = N/ (x(7o) + 1) we can write the Kovchegov equation 

as 

^n+l = (X(70) + - ^n). (7) 

One can recognize ((Tj) as the famous logistic map jl^; u = x(7o) + 1 is called the Malthusian 
parameter. The properties of its solutions are listed below (see e. g. .16]). 

1. 1 < c<j < 3: the scattering amplitude folows the sigmoid curve towards the saddle point 
value Dl/ (the fixed point). Since D^/ < 1, the black disc limit is not achieved even at 
asymptotically high energies. 

2. 3 < uj < UJ2, where 002 ~ 3.449. At these values of u the scattering amplitude at 
n — s> 00 does not converge to a single limit - instead, it oscillates between two fixed 
points. Mathematically, 91 develops a pitchfork bifurcation at = 3 at which point 
the 2-cycle begins. As seen in Fig. |21the next bifurcation happens at UJ4 ~ 3.449; it 
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FIG. 2: Bifurcation map: the fixed points of small x evolution as a function of Malthusian param- 
eter (J. 

starts a 4-cycle. The 8-cycle starts at cjg ^ 3.56 and so on. Note, that the behavior of 
91 at n — > oo is independent of the inital condition at n = 1 as long as u; < 3.57. 

3. liu> ciJap, where u!a_p ~ 3.57 is the accumulation point, periodicity gives way to chaos. 
In other words the scattering amplitude exhibits irregular, unpredictable behavior 
which manifests itself in a sensitivity to small changes in the initial condition. Recall 
that the BFKL saddle point is 7bfkl = 1/2 at which x(7bfkl) = 4 In 2 ^ 2.77, thus 
i^BFKL = 3.77. In perturbation theory u cannot take values less than cubfkl- Therefore, 
we conclude that the perturbative high energy evolution is chaotic. 

By averaging over all events one can define the mean value of the scattering amplitude. 
However, this procedure hides a lot of interesting physics. The most obvious example of this 
is diffraction, which measures the strength of fluctuations in the inelastic cross section jl7|. 
The arguments given in points 2-3 above show that diffraction is a significant part of the 
total inelastic cross section at very high energies, and is universal (independent from the 
properties of the target). 

It is important to emphasize, that our model treats the high energy evolution process 
classically. We neglected the fact that the gluon emission is a stochastic process. The 
emission "time", i.e. the rapidity interval dy over which a gluon is radiated, varies from 
event to event. In other words, gluon emission is a quantum process which may or may not 
occur with a certain probability once the rapidity interval of the collision y is increased by 
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dy. Full treatment of the discrete BK equation requires taking these effects into account. 
However, unfortunately BK equation is known to resist all attempts of analytical solution, 
and our hope at present is to develop a meaningful approximation. Thus, in our paper we 
suggested an approximation in which the gluons are emitted over a fixed "time" defined 
by agdy = C with C = 1. To justify this assumption, let us note that BFKL takes into 
account only fast gluons, i.e. those with C ~ 1. It is beyond the leading logarithmic (LL) 
approximation to take into account slow gluons. Moreover, it^ known that an account of 
NLL corrections effectively leads to imposing a rapidity veto [18] on the emission of gluons 
with close rapidities, which restricts production of gluons with small C (this is due to an 
effective repulsion between the emitted gluons induced at the NLL level). Therefore C is 
bounded from below by a number close to one. On the other hand the probability that no 
gluon is emitted when C becomes larger than one is very small if we choose the high density 
initial condition, such as the one given by the McLerran-Venugopalan modelj^. Therefore, 
C takes random values around 1, but the effective dispersion can be expected quite small. 

We realize that the quantum fluctuations can affect the effective value of C and push the 
onset of chaos to a different kinematic region, but we believe that this effect is not going to 
be dramatic. The present paper is only the first step in the exploration of the discrete BK 
equation. Our aim is to stress its highly nontrivial structure which might have an important 
impact on the high energy theory and phenomenology. The question about the stability of 
the found peculiar classical solution with respect to the stochastic quantum fluctuations is 
very important and we are going to address it in the forthcoming work. 

The model used in this letter is admittedly oversimplified: we neglected the diffusion 
in transverse momentum, stochasticity of gluon emission and the dynamical fluctuations 
beyond the mean field approximation. Nevertheless, we hope that at least some of the 
features of discrete quantum evolution at small x will survive a more realistic treatment. 
The chaotic features of small x evolution open a new intriguing prospective on the studies 
of hadron and nuclear interactions at high energies. 

Acknowledgments 

The authors would like to thank Yuri Kovchegov, Alex Kovner, Misha Kozlov, Genya 
Levin, Larry McLerran, Al Mueller, Anna Stasto for informative and helpful discussions 



6 



and comments. This research was supported by the U.S. Department of Energy under 
Grant No. DE-AC02-98CH10886. 



[1] L. V. Gribov, E. M. Levin and M. G. Rysliin, Phys. Rept. 100, 1 (1983). 
[2] A. H. Mueller and J. w. Qiu, Nucl. Phys. B 268, 427 (1986). 

[3] L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 3352 (1994) |arXiv:hep-ph/9311205| ; 



Phys. Rev. D 49, 2233 (1994) |arXiv:hep-ph/9309289, ; Phys. Rev. D 50, 2225 (1994) 
|arXiv:hep-ph/9402335] . 

[4] J. Jalihan-Marian, A. Kovner, A. Leonidov and H. Weigert, Phys. Rev. D 59, 014014 (1999) 
|arXiv:hep-ph/9 706377 . J. Jalilian-Marian, A. Kovner and H. Weigert, Phys. Rev. D 59, 
014015 (1999) arXiv:hep-ph/9709432l. 

[5] E. lancu, A. Leonidov and L. D. McLerran, Nucl. Phys. A 692, 583 (2001) 
|arXiv:hep-ph/0011241| ; Phys. Lett. B 510, 133 (2001) |arXiv:hep-ph/0102009| . E. lancu and 



L. D. McLerran, Phys. Lett. B 510, 145 (2001) ar Xiv:hep-ph /0103032 . 
[6] Y. V. Kovchegov, Phys. Rev. D 54, 5463 (1996) |arXiv:hep-ph/9605446 ; Phys. Rev. D 55, 
5445 (1997) arXiv:hep-ph/9701229 . 



[7] E. lancu, A. H. Mueher and S. Munier, arXiv:hep-ph/0410018 



[8] E. lancu and A. H. Mueher, Nucl. Phys. A 730, 494 (2004) |arXiv:hep-ph/0309276| . 
[9] E. lancu and D. N. Triantafyllopoulos, arXiv:hep-ph/0411405l 



[10] A. H. Mueher, A. I. Shoshi and S. M. H. Wong, arXiv:hep-ph/0501088 



[11] I. Bahtsky Nucl. Phys. B 463, 99 (1996) |arXiv:hep-ph/9 509348] . 

[12] Y. V. Kovchegov, Phys. Rev. D 60, 034008 (1999) |arXiv:hep-ph/9901281| . 

[13] E. A. Kuraev, L. N. Lipatov and V. S. Fadin, Sov. Phys. JETP 45, 199 (1977) [Zh. Eksp. 

Teor. Fiz. 72, 377 (1977)]; 1. I. Bahtsky and L. N. Lipatov, Sov. J. Nucl. Phys. 28, 822 (1978) 

[Yad. Fiz. 28, 1597 (1978)]. 
[14] Y. V. Kovchegov, Phys. Rev. D 61, 074018 (2000) |arXiv:hep-ph/990521 4|. 
[15] J. von Neumann, Applied Math. Series 12, National Bureau of Standards, 1951, 3638. 
[16] A concise review on the properties of logistic map and the logistic equation can be found at 

|http:/ / mathworId.wolfram.com/LogisticMap.html 



[17] E.L. Feinberg and I.Ya. Pomeranchuk Suppl. Nuovo Cim. 3 (1956), p. 652; 



7 



M. L. Good and W. D. Walker, Phys. Rev. 120, 1857 (1960); 
J. Pumplin Phys. Rev. D8, 2899 (1973). 
[18] C. R. Schmidt, Phys. Rev. D 60, 074003 (1999) |arXiv:hep-ph/9901397| ; 

J. R. Forshaw, D. A. Ross and A. Sabio Vera, Phys. Lett. B 455, 273 (1999) 
|arXiv:hep-ph/9903390] ; 

G. Chachamis, M. Lubhnsky and A. Sabio Vera, Nucl. Phys. A 748, 649 (2005) 
|arXiv:hep^h/0408333 . 



